Millimeter-wave  Interferometric  Sensing 


by  Charles  R.  Dietlein,  Amir  Zaghloul,  and  David  A.  Wikner 


ARL-TR-5668 


September  2011 


Approved  for  public  release;  distribution  unlimited. 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position 
unless  so  designated  by  other  authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an  official  endorsement  or 
approval  of  the  use  thereof. 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


Army  Research  Laboratory 

Adelphi,MD  20783-1 197 


ARL-TR-5668  September  2011 


Millimeter-wave  Interferometric  Sensing 


Charles  R.  Dietlein,  Amir  Zaghloul,  and  David  A.  Wikner 
Sensors  and  Electron  Devices  Directorate,  ARL 


Approved  for  public  release;  distribution  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the 
burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

This  report  explores  the  feasibility  of  short-range  ground-based  millimeter-wave  sensing  systems  based  on  interferometric 
imaging  techniques.  In  comparison  with  infrared  and  visible  wavelength  focal-plane  arrays  (FPA),  millimeter-wave  FPAs 
always  have  fewer  detectors  and  image  pixels,  poorer  spatial  resolution,  lower  signal-to-noise  ratio,  and  an  increase  of 
illumination-based  artifacts  such  as  glint  in  the  images.  Calibration  and  scanning  schemes  induce  additional  complications. 

These  difficulties  are  ultimately  a  consequence  of  the  millimeter-wave’s  much  longer  wavelengths  compared  to  visible  or 
infrared  frequencies.  This  research  effort  addresses  a  specific  limitation:  low-detector-count  FPAs.  An  alternative  approach  to 
traditional,  high-detector-count,  staring  FPAs  is  interferometric  (coherence)  imaging,  wherein  an  array  with  a  low  number  of 
coherent  detectors  is  exploited  to  reconstruct  the  field  distribution  in  the  object  plane.  This  interferometric  approach  is  used 
extensively  by  radio  astronomers  at  lower  frequencies.  In  this  report,  we  study  this  technique's  validity  at  220  GFlz  for 
imaging  terrestrial  scenes. 

15.  SUBJECT  TERMS 

Millimeter-wave,  imaging,  coherence,  interferometry,  standoff 

17.  LIMITATION  18.  NUMBER 
OF  OF 

ABSTRACT  PAGES 

UU  20 

Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


19a.  NAME  OF  RESPONSIBLE  PERSON 

Charles  R.  Dietlein 

19b.  TELEPHONE  NUMBER  ( Include  area  code ) 

(301)  394-1523 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

b.  ABSTRACT 

C.  THIS  PAGE 

Unclassified 

Unclassified 

Unclassified 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

September  2011  Final 

4.  TITLE  AND  SUBTITLE 

Millimeter-wave  Interferometric  Sensing 


6.  AUTHOR(S) 

Charles  R.  Dietlein,  Amir  Zaghloul,  and  David  A.  Wikner 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
ATTN:  RDRL-SER-M 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED  (From  -  To) 

March  to  June  2010 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-TR-5668 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


11 


Contents 


List  of  Figures  iv 

Summary  1 

1.  Introduction  3 

2.  Methods,  Assumptions,  and  Procedures  3 

3.  Results  and  Discussion  6 

4.  Conclusions  11 

5.  References  12 

Distribution  List  13 


iii 


List  of  Figures 


Figure  1 .  Geometry  illustrating  the  van  Cittert-Zernike  theorem  discussed  in  the  text . 4 

Figure  2.  Plot  of  phase  difference  between  detectors  for  various  array  sizes  (shown  in  the 
legend  in  units  of  meters),  as  a  function  of  range  between  the  array  and  object  planes. 
Frequency  was  fixed  at  220  GHz.  The  phase  differences  <j>g  «  1  imply  a  paraxial 
approximation.  The  black  square  indicates  where  the  experiments  and  simulations  in 
this  report  were  perfonned . 6 

Figure  3.  (Left)  Array  geometry  showing  32  detector  positions  arranged  in  a  spiral 

configuration,  and  (right)  sampling  in  the  uv  plane . 7 

Figure  4.  (Left)  A  histogram  containing  the  uv  plane  coverage  shown  in  Cartesian 

coordinates  in  figure  3  (right).  The  histogram  contains  496  samples  in  50  bins.  (Right)  A 
histogram  showing  uv  plane  coverage  for  a  rectangularly-gridded  array,  with 
approximately  the  same  number  of  detectors  and  physical  size  as  that  of  figure  3  (left) . 8 

Figure  5.  (Left)  Reconstruction  of  the  simulated  point  source  at  a  range  of  0.5  m  from  the 
spiral  array,  and  (right)  best-fit  2-D  Gaussian  (opaque  surface)  with  the  reconstructed 
object  plane  (mesh  surface),  for  purposes  of  determining  the  FWHM  of  the  point  source . 9 

Figure  6.  (Left)  Reconstruction  of  the  two  simulated  point  sources  at  a  range  of  0.5  m  from 
the  spiral  array,  and  (right)  best-fit  2-D  Gaussians  (opaque  surface)  with  the 
reconstructed  object  plane  (mesh  surface) . 9 

Figure  7.  Measured  electric  field  amplitude  pattern  of  two  waveguide  aperture  sources. 

Range  to  the  sources  is  approximately  0.5  m,  and  they  are  separated  by  approximately 
6  cm.  The  fringes  are  tilted  slightly,  as  the  sources  were  not  perfectly  aligned  with  the  y- 
axis.  The  nominal  positions  of  the  modeled  detectors  are  depicted  by  the  overlaid  dots . 10 

Figure  8.  (Left)  Measured  data,  reconstructed  electric  field  distribution  at  object  plane  and 

(right)  simulation  results  with  geometry  matching  the  experimental  setup . 1 1 


IV 


Summary 


Millimeter-wave  sensing  systems  designed  for  imaging  almost  always  produce  inferior  images 
compared  to  their  infrared  and  visible  wavelength  counterparts  for  a  variety  of  fundamental 
reasons.  These  reasons  include  low  detector/pixel  count,  poor  signal-to-noise  ratio,  artifacts 
introduced  by  illumination,  and  issues  arising  from  calibration  or  scanning  schemes.  This 
research  effort  is  intended  to  address  the  issue  of  low-detector-count  focal-plane  arrays  (FPA). 
Radio  astronomers  have  dealt  with  this  issue  for  many  decades  by  using  interferometry,  or 
coherence  imaging.  This  report  details  how  interferometric  imaging  techniques  can  be  used  in 
the  millimeter-wave  regime  for  terrestrial  applications. 
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1.  Introduction 


This  research  effort  explores,  through  computer  simulations  and  experimental  verification,  the 
viability  of  interferometric  techniques  for  millimeter-wave  imaging.  A  fundamental  limitation  in 
millimeter- wave/terahertz  imaging  is  the  lack  of  high-detector-count  FPAs.  However,  in  the 
radio  frequency  through  terahertz  regime,  coherent  (phase-sensitive)  detectors  are  available. 
Infrared  and  visible  images  are  acquired  with  direct  (amplitude-only)  detectors.  Coherent 
measurement  of  the  incident  electric  field  provides  additional  infonnation  regarding  the  spatial 
amplitude  distribution  in  the  object  plane.  The  scope  of  this  research  is  to  explore  an 
interferometric  imaging  array  geometry  with  the  goal  of  demonstrating  that  a  sparsely  populated 
coherent  array  can  achieve  comparable  spatial  resolution  to  that  achieved  by  a  traditional  direct- 
detection  fully  populated  FPA  with  similar  physical  dimensions.  The  computer  program 
MATLAB®  is  employed  to  simulate  the  image  amplitude  distribution  and  quantify  the  spatial 
resolution  of  the  array,  and  a  table-top  experiment  is  configured  to  verify  the  simulations. 


2.  Methods,  Assumptions,  and  Procedures 


An  imaging  interferometer  consists  of  N  coherent  detectors  arranged  in  an  array;  we  will 
consider  planar  arrays  for  simplicity.  Each  detector  measures  the  amplitude  and  phase  of  the 
incident  radiation,  while  each  pair  of  detectors  measures  one  spatial  Fourier  component  of  the 
incident  radiation.  The  specific  spatial  component  is  defined  by  the  pair’s  physical  separation, 
called  a  baseline.  For  N  detectors,  there  are 


P  = 


v2y 


(1) 


unique  2-combinations.  Each  spatial  Fourier  component  is  represented  by  a  single  point  in  the 
Fourier  plane,  known  as  the  uv  plane.  Image  quality  is  usually  determined  by  the  uv  plane 
sampling.  Since  each  baseline  measures  a  single  spatial  frequency,  the  optimum  array  geometry 
contains  no  redundant  spacings.  Typical  geometries  include  a  logarithmically  spaced  one¬ 
dimensional  (1-D)  array,  or  a  spiral  two-dimensional  (2-D)  array.  While  image  quality  is 
primarily  dependent  on  uv  plane  sampling,  it  can  be  improved  via  standard  reconstruction  post¬ 
processing  algorithms  that  are  used  by  radio  astronomers. 

The  van  Cittert-Zemike  theorem  explains  the  fundamentals  of  interferometric  imaging.  A  more 
detailed  version  of  this  discussion  can  be  found  in  many  textbooks  (1 ).  Briefly,  a  source  of 
electromagnetic  radiation  occupies  the  infinitesimal  surface  dS'  on  the  z  =  0  plane,  at  xy 
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coordinates  x  =  x'  and  y  =  y'.  Aq  is  the  amplitude  of  the  source.  An  array  of  detectors  is  in  the  xy 
plane  at  z  =  Zq.  This  geometry  is  illustrated  in  figure  1 . 


For  a  pair  of  detectors  at  (xi,  y\,  Zq)  and  (X2,  yi,  Zq),  the  contributions  to  the  total  electric  field 
provided  by  the  source  d5"  can  be  expressed  by 


and 


dZT,  =  (A0  /  rt  )exp  {jot  -  jkrx  )dS' 


dE2  =  (A0  /  r2  )exp  ( jcot  -  jkr2  )dS" 


(2) 


(3) 


where  r.  =  [(x'-xy )2  +  [y'-y j  )2  +  Z, 


1/2 


and  k  =  2n  I Z.  The  mutual  coherence  function  of  the 


electric  field  at  detector  positions  (xi,  y i,  Zo)  and  (x2,  yi,  Zq)  is  given  by 


Q,=J 


r7/,(x',y')cxp[/7c(/| 


(4) 


where  cr£  is  the  object  plane  amplitude  distribution,  and  in  the  paraxial  approximation, 

(xf  +  vf)-(x22  +yl)  (x,  -  x2  )x'+(yx  -  y2 )/ 


27  7 


(5) 


The  phase  difference  between  a  pair  of  detectors  in  a  planar  array,  due  to  wave-front  curvature,  is 
given  by 


$8  — 


k\x{  +  y\ )-  (x2  +  y2  )] 


2Zn 


(6) 
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If  we  define  the  baselines  u  =  xx  -  x2  and  v  =  yl  -  y2 ,  along  with  the  angular  coordinates 
£  =  jc'/Z0  ,  T)  =  y  /Z0  and  the  approximation  1/  r,  r2  «  Z02 ,  the  mutual  coherence  function  can  be 
written  as  (2) 

oo  oo 

C\t2  (U’  V)  =  exp  (/ (f)s  )  J  |  (£,//)  exp  [-  jk  (u  g  +  v//  )]d  cd  // .  (7) 

-oo-oo 

When  the  phase  shift  exp (j<f>e)  can  be  ignored,  the  mutual  coherence  function  above  relates  the 

coherence  in  the  array  plane  to  the  amplitude  distribution  in  the  object  plane  via  a  Fourier 
transform.  When  the  phase  shift  is  significant,  the  correction  based  on  the  array  geometry  and 
distance  to  the  object  plane  can  be  performed  by  adding  a  tenn  to  the  measured  phase  at  each 
detector,  before  computing  the  correlation  function  per  baseline,  C,  =  A ,  exp  (jA<f>l ) .  In  the 

correlation  function,  Al  =  EmEn  is  the  product  of  the  measured  field  amplitudes  at  detectors  m 
and  n  ,  and  A0,  =  (f>m  -(f) n  is  the  difference  in  phase  between  the  two  detectors. 

For  common  array  sizes  in  the  millimeter- wave/terahertz  frequency  range,  imaging  at  realistic 
ranges,  the  far-field  condition  of  <j>s  « 1  is  not  met,  and  <f)s  must  be  adjusted  in  software  after 

measurement.  The  relationship  between  range,  array  size,  and  phase  shift  is  shown  in  figure  2. 
One  may  ask  why  we  use  the  van  Cittert-Zemike  theorem,  strictly  valid  only  in  the  paraxial 
approximation,  when  we  are  in  at  least  one  of  the  three  near-field  conditions  (3).  The  reason  is 
that  the  paraxial  approximation  provides  a  useful  form  of  the  coherence  function,  because  aE 

can  be  reconstructed  easily  through  Fourier  inversion.  Reconstruction  of  the  amplitude 
distribution  from  the  correlation  function  of  each  baseline  is  perfonned  by 

oE  (£  7i)=Yj  MC/  )cos  [k(u£  +  vy/)]- ~s{C/  )sin  [k{u,4  +  vy/)]} . 

/=! 
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Figure  2.  Plot  of  phase  difference  between  detectors  for  various  array  sizes  (shown 

in  the  legend  in  units  of  meters),  as  a  function  of  range  between  the  array  and 
object  planes.  Frequency  was  fixed  at  220  GHz.  The  phase  differences  <j>s  «  1 
imply  a  paraxial  approximation.  The  black  square  indicates  where  the 
experiments  and  simulations  in  this  report  were  performed. 


3.  Results  and  Discussion 


We  focused  our  efforts  around/=  220  GHz,  primarily  because  it  is  one  of  the  few  bands  with  an 
acceptable  atmospheric  transmission  window  above /=  100  GHz.  Additionally,  the  wavelength 
of  approximately  X  =  1 .4  mm  provides  adequate  spatial  resolution  for  a  variety  of  tasks  with 
moderate  array  sizes.  Finally,  it  is  easy  to  find  commercial  hardware  that  operates  at  this 
frequency.  While  this  approach  will  work  at  arbitrary  ranges,  we  constrained  our  simulations  to 
those  that  could  be  replicated  on  a  laboratory  bench  with  off-the-shelf  equipment.  We  chose  a 
range  of  approximately  0.5  m  and  an  array  size  of  approximately  20  mm  .  Figure  3  (left)  shows 
one  of  the  simulated  array  geometries  we  explored,  a  planar  Archimedean  spiral,  with  N  =  32 
detectors.  The  array’s  sampling  in  the  uv  plane  is  shown  in  figure  3  (right). 
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Figure  3.  (Left)  Array  geometry  showing  32  detector  positions  arranged  in  a  spiral  configuration,  and  (right) 
sampling  in  the  uv  plane. 

While  figure  3  (right)  appears  to  be  sampling  the  Fourier  space  in  a  non-uniform  way,  it  is 
difficult  to  visually  determine  if  there  is  redundant  sampling  at  similar  spatial  frequencies.  A 
simple  check  can  be  provided  via  a  histogram  of  the  baseline  lengths,  which  appears  in  figure  4 
(left).  There  are  p  =  496  baselines  in  this  array,  and  the  histogram  contains  50  bins.  It  can  be  seen 
that  this  geometry  is  neither  a  poor  choice  (e.g.,  a  rectangular  gridded  array,  which  samples 
relatively  few  spatial  frequencies)  nor  an  ideal  one  (in  which  the  histogram  would  be  nominally 
flat  from  the  minimum  spatial  frequency  to  the  maximum).  There  are  exactly  490  unique 
baseline  lengths  in  the  array;  a  non-redundancy  figure  of  merit  (FOM)  can  be  defined  as  the 
number  of  unique  baseline  lengths  divided  by  the  total  number  of  baselines  in  the  array.  For  this 
array,  that  FOM  is  0.988.  For  comparison,  figure  4  (right)  contains  a  histogram  of  the  baseline 
lengths  for  a  rectangular- gridded  array  with  N  =  36  detectors  (and  therefore  p  =  630  baselines) 
and  the  same  approximate  physical  size  as  the  planar  spiral  array;  inset  is  the  rectangular  array 
configuration  ( cf.  figure  3  (left)).  Straightforward  geometry  shows  that  there  are  exactly  19 
unique  baseline  lengths  in  the  rectangular-gridded  array.  The  non-redundancy  FOM  for  this  array 
is  0.030. 
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Figure  4.  (Left)  A  histogram  containing  the  uv  plane  coverage  shown  in  Cartesian  coordinates  in  figure  3  (right). 
The  histogram  contains  496  samples  in  50  bins.  (Right)  A  histogram  showing  uv  plane  coverage  for  a 
rectangularly-gridded  array,  with  approximately  the  same  number  of  detectors  and  physical  size  as  that 
of  figure  3  (left). 

The  spatial  resolution  of  an  interferometric  array  can  be  approximated  by  XZ^Ib,  where  b  is  the 
largest  baseline  in  the  array.  For  the  array  shown  in  figure  3  (left),  b  ~  19  mm.  At  a  range  of 
0.5  m  the  spatial  resolution  should  be  approximately  3.6  cm.  Figure  5  (left)  shows  a 
reconstructed  amplitude  distribution  of  a  point  source  at  0.5  m.  To  verify  the  full-width  half- 
maximum  (FWHM)  spot  size,  a  fitting  code  using  the  Nelder-Mead  direct  search  simplex 
algorithm  ( 4 )  was  written  to  find  the  best-fit  2-D  Gaussian  to  the  data.  The  general  form  of  the 
2-D  Gaussian  is  used:  /[x,  y)  =  A  exp-[u(x  -  xo)~  +  2 b(x  -  xo)(y  -yo)  +  c(y  -yo)  ],  where  a,  b,  and 
c  are  functions  of  both  the  rotation  angle,  and  x  and  y  Gaussian  widths.  Figure  5  (right)  shows  the 
results  of  the  fitting  algorithm.  The  mesh  surface  is  the  simulated  data,  while  the  opaque  surface 
is  the  best-fit  2-D  Gaussian.  The  mean  of  the  x  and  y  Gaussian  widths  was  a  ~  0.03  m; 
multiplication  by  x/21n  2  results  in  a  FWHM  of  3.55  cm,  which  is  very  close  to  the  expected 
3.6  cm. 
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Figure  5.  (Left)  Reconstruction  of  the  simulated  point  source  at  a  range  of  0.5  m  from  the  spiral  array,  and 

(right)  best-fit  2-D  Gaussian  (opaque  surface)  with  the  reconstructed  object  plane  (mesh  surface),  for 
purposes  of  determining  the  FWHM  of  the  point  source. 

Next,  we  investigated  multiple  point  sources.  With  the  same  planar  Archimedean  spiral  array  and 
distance  to  the  object  plane,  we  added  a  second  point  source  6  cm  above  the  first.  The  results  in 
figure  6  show  that  the  two  sources  are  clearly  resolved.  A  similar  best- fit  Gaussian  code  was 
developed  to  accommodate  fitting  to  multiple  Gaussians  and  was  applied  to  the  two-source 
simulation  results. 


Figure  6.  (Left)  Reconstruction  of  the  two  simulated  point  sources  at  a  range  of  0.5  m  from  the  spiral  array,  and 
(right)  best-fit  2-D  Gaussians  (opaque  surface)  with  the  reconstructed  object  plane  (mesh  surface). 
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To  determine  if  the  simulation  results  were  valid,  we  constructed  a  simple  table-top  experiment 
consisting  of  a  220-GHz  transmitter  and  receiver.  We  installed  a  waveguide  directional  coupler 
at  the  transmitter  output,  and  with  a  pair  of  attenuators,  set  the  output  power  at  both  waveguides 
(coupled  and  through)  to  be  as  close  to  identical  as  possible.  This  measurement  was  performed 
using  a  calibrated  substitution  power  meter.  The  spacing  between  the  waveguide  apertures  was 
approximately  6  cm.  We  raster-scanned  the  receiver  over  the  spatial  region  where  the  32 
detectors  were  modeled.  The  field  was  sampled  on  a  square  grid  with  a  pixel  dimension  of 
0.794  tnm  per  side.  The  results  are  shown  in  figure  7  (amplitude  only,  phase  not  shown),  with  the 
nominal  detector  positions  of  the  modeled  Archimedean  spiral  overlaid.  The  interference  pattern 
is  slightly  tilted  because  the  waveguide  apertures  were  not  perfectly  aligned  with  the  y-axis  of 
the  experiment. 


Figure  7.  Measured  electric  field  amplitude  pattern  of  two  waveguide  aperture  sources. 

Range  to  the  sources  is  approximately  0.5  m,  and  they  are  separated  by  approximately 
6  cm.  The  fringes  are  tilted  slightly,  as  the  sources  were  not  perfectly  aligned  with 
the  y-axis.  The  nominal  positions  of  the  modeled  detectors  are  depicted  by  the 
overlaid  dots. 

The  measured  complex  electric  field  was  interpolated  to  a  grid  with  three  times  more  points  in  x 
andy,  and  then  the  electric  field  amplitude  and  phase  were  interpolated  to  the  nearest  points  of 
the  modeled  spiral.  This  produces  a  1x32  vector  with  a  measured  electric  field  amplitude  and 
phase  at  each  detector’s  position  in  the  modeled  spiral.  This  1x32  complex  vector  was  loaded 
into  the  MATLAB  code  used  in  the  previous  simulations,  and  the  code  was  started  at  the 
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subroutines  that  process  the  electric  field  amplitude  and  phase  from  each  detector.  The  next  step 
in  the  processing  is  to  calculate  the  correlation  function  for  each  baseline  in  the  spiral  as 
mentioned  earlier,  C,  =  A,  exp  (jA(f)l ) ,  where  A,  =  EmEn  is  the  product  of  the  measured  field 

amplitudes  at  detector  locations  in  and  n  ,  and  Atf>l  =  (f>m  -  (f>n  is  the  difference  in  phase  between 
the  two  detector  locations.  The  reconstruction  result  is  shown  in  figure  8  (left). 
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Figure  8.  (Left)  Measured  data,  reconstructed  electric  field  distribution  at  object  plane  and  (right) 
simulation  results  with  geometry  matching  the  experimental  setup. 

Figure  8  (right)  contains  the  results  of  the  same  geometry,  but  with  purely  simulated  data.  The 
agreement  between  simulation  and  experiment  is  very  good,  and  verifies  that  the  concepts 
explained  in  section  3  are  valid. 


4.  Conclusions 


We  explored  and  verified  the  feasibility  of  interferometric  imaging  techniques  for  use  at 
millimeter- wave  frequencies  and  short  ranges.  The  paraxial  approximation  of  the  van  Cittert- 
Zernike  theorem  that  radio  astronomers  use  is  applicable  in  the  situations  described  by  this 
report,  as  long  as  the  phase  front  curvature  is  compensated  for.  This  is  equivalent  to  re-focusing 
the  array  to  the  near-field,  a  technique  that  has  been  discussed  by  some  in  the  radio  astronomy 
community  (5)  for  studying  astronomical  objects  in  our  solar  system. 
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